Effect of the particle-hole channel on BCS-Bose-Einstein condensation crossover in atomic Fermi 
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BCS-Bose-Einstein condensation (BEC) crossover is effected by increasing pairing strength between 
fermions from weak to strong in the particle-particle channel. Here we study the effect of the particle-hole 
channel on the zero T gap A(0), superfluid transition temperature Tc and the pseudogap at Tc, as well as the 
mean- field ratio 2A(0) /T^^, from BCS through BEC regimes, in the framework of a pairing fluctuation theory 
which includes self-consistently the contributions of finite-momentum pairs. These pairs necessarily lead to a 
pseudogap in single particle excitation spectrum above and below Tc. We sum over the infinite particle-hole 
ladder diagrams so that the particle-particle and particle-hole T-matrices are entangled with each other. We find 
that the particle-hole susceptibility has a complex dynamical structure, with strong momentum and frequency 
dependencies, and is sensitive to temperature, gap size and interaction strength. We conclude that neglecting 
the self-energy feedback causes a serious over-estimate of the particle-hole susceptibility. In the BCS limit, 
the particle-hole channel effect may be approximated by the same reduction in the overall pairing strength so 
that the ratio 2A(0)/rc is unaffected, in agreement with Gor'kov et al. to the leading order. However, the 
effect becomes more complex and pronounced in the crossover regime, where the particle-hole susceptibility 
is reduced by both a smaller Fermi surface and a big (pseudo)gap. Deep in the BEC regime, the particle-hole 
channel contributions drop to zero. We propose that precision measurements of the magnetic field for Feshbach 
resonance at low temperatures as a function of density can be used to quantify the particle-hole susceptibility 
and test different theories. 

PACS numbers: 03.75. Ss, 03.75.Nt, 74.20.-z, 74.25.Dw 



I. INTRODUCTION 

BCS-Bose-Einstein condensation (BEC) crossover has 
been an interesting research subject since 1980's [1-29]. The 
experimental realization of BCS-BEC crossover in ultracold 
atomic Fermi gases [30-39], with the help of Feshbach reso- 
nances, has given it a strong boost [40-59] over the past sev- 
eral years. When the pairing interaction is tuned from weak 
to strong in a two component Fermi gas, the superfluid behav- 
ior evolves continuously from the type of BCS to that of BEC 
[2, 3, 29]. 

In such a fundamentally fermionic system, superfluidity 
mainly concerns pairing, namely, interactions in the particle- 
particle channel. In contrast, the particle-hole channel mainly 
causes a chemical potential shift, and is often neglected [60]. 
For example, in a conventional superconductor, the chemical 
potential below and above Tc are essentially the same, and 
thus its dependence on the temperature and the interaction 
strength has been completely neglected in the weak coupling 
BCS theory for normal metal superconductors. On the other 
hand, Gor'kov and Melik-Barkhudarov (GMB) [61] consid- 
ered the lowest order correction from the particle-hole chan- 
nel, (which has been referred to as induced interaction in the 
literature), and found that both Tc and zero temperature gap 
A(0) are suppressed by a big factor of (4e)^/'^ w 2.22. Berk 
and Schrieffer [62] also studied a similar effect in the form 
of ferromagnetic spin correlations in superconductors. De- 
spite the big size of the GMB correction, the effect of the 
particle-hole channel has been largely overlooked in the theo- 
retical study of BCS-BEC crossover, until it has become real- 
istic to achieve such crossover experimentally in atomic Fermi 



gases. Heiselberg and coworkers [63] considered the effect 
of the lowest order induced interaction in dilute Fermi gases 
and generalized it to the case of multispecies of fermions as 
well as the possibility of exchange of bosons. Kim et al. 
[64] considered the lowest order induced interactions in op- 
tical lattices. Within the mean-field treatment and without in- 
cluding the excitation gap in the particle and hole propaga- 
tors, these authors found the same effective overall interaction 
at zero T and at and hence an unaffected mean-field ratio 
2A(0)/fcBT'c « 3.53. Martikainen ef fl/. [65] considered the 
lowest order induced interactions in a three-component Fermi 
gas. It has become clear that including only the perturbative 
lowest order induced interaction is not appropriate away from 
the weak coupling BCS regime. Yin and coworkers [66] went 
beyond the lowest order and considered the induced interac- 
tions from all particIe-hoIe ladder diagrams, i.e., the entire 
particle-hole T-matrix. However, in all the above works, only 
the bare particIe-hoIe susceptibility Xph was considered, and 
it was averaged on-shell and only on the Fermi surface, with 
equal momenta for the particle and the hole propagators. No 
self-energy feedback was included. Therefore, there was nec- 
essarily no pseudogap in the fermion excitation spectrum at 
Tc. This is basically equivalent to replacing the particle-hole 
susceptibility Xph essentially temperature independent 

constant, leading to a simple downshift in the pairing interac- 
tion. 

As the gap and Tc increase with interaction strength, it can 
naturally be expected that the contribution from the particle- 
hole channel, or the induced interaction, will acquire a signifi- 
cant temperature dependence. More importantly, the presence 
of the (pseudo)gap serves to suppress the particle-hole fluctu- 
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ations (which tend to break pairing). In other words, neglect- 
ing the feedback of the gap related self energy in the particle- 
hole susceptibility is expected to cause an over-estimate of the 
particle-hole channel contributions. Therefore, a proper treat- 
ment should include the gap effect in the particle-hole suscep- 
tibility. In addition, the lowest order treatment is no longer 
appropriate away from the weak coupling regime. 

Furthermore, it has now been established that as the pairing 
interaction increases, pseudogap develops naturally [29, 67]. 
Experimental evidence for its existence comes from high Tc 
superconductors [29, 68-70] as well as atomic Fermi gases 
[71-75]. Therefore, a theory with proper treatment of the 
pseudogap effect is necessary in order to arrive at results that 
can be quantitatively compared with experiment. For the same 
reason, the effect of the particle-hole channel needs also to be 
studied within such a theory. 

In this paper, we explore the particle-hole channel effect 
based on a pairing fluctuation theory [21, 22], originally de- 
veloped for treating the pseudogap phenomena of high Tc su- 
perconductors. This theory has been successfully applied to 
atomic Fermi gases and has been generating results that are 
in good agreement with experiment [29, 71, 73]. Here we 
include the entire particle-hole T-matrix, with gap effect in- 
cluded in the fermion Green's functions. Instead of a simple 
average of the particle-hole susceptibility Xph on the Fermi 
surface, here we choose to average at two different levels - 
one on the Fermi surface, one over a narrow momentum shell 
around the Fermi level. We find that Xph has very strong fre- 
quency and momentum as well as temperature dependencies. 
It is sensitive to the gap size. Therefore, self-consistently in- 
cluding the self-energy feedback is important. For both levels 
of average, we find that while in the BCS Umit, the particle- 
hole channel effect may be approximated by a downshift in the 
pairing strength so that the ratio 2A(0)/rc is unaffected, the 
situation becomes more complex as the interaction becomes 
stronger where the gap is no longer very small. Significant dif- 
ference exists for these two levels of averaging. The particle- 
hole susceptibility is reduced by both a smaller Fermi surface 
and a big (pseudo)gap in the crossover regime. Deep in the 
BEC regime, the particle-hole channel contributions drop to 
zero. Without including the incoherent part of the self energy, 
we find that at unitarity, TcZ-Bp « 0.217, in reasonable agree- 
ment with experiment. 

The rest of this paper is arranged as follows. In Sec. II, 
we first give a brief summary of the pairing fluctuation theory 
without the particle-hole channel effect. Then we derive the 
theory with particle-hole channel included, starting by study- 
ing the dynamic structure of the particle-hole susceptibility. 
Next, in Sec. Ill, we present numerical results, showing the 
effect of the particle-hole channel on the zero T gap, transition 
temperature Tc and pseudogap at Tc, as well as the mean-field 
ratio 2A(0)/Tj?^^. We also discuss and compare our value of 
Tc/Ep with experiment and those in the literature. Finally, 
we conclude in Sec. IV. More detailed results on the dynamic 
structure of the particle-hole susceptibility are presented in the 
Appendix. 



] + L 



+ 



S = 




" ~ V ^sc 



Figure 1. Feynman diagrams for the particle-particle channel T- 
matrix tpg and the self energy 'E{K). The dotted lines represent 
the bare pairing interaction U. The dashed line, tx, represents the 
superfluid condensate. 



II. PAIRING FLUCTUATION THEORY WITH THE 
PARTICLE-HOLE CHANNEL EFFECT INCLUDED 

A. Summary of the pairing fluctuation theory without the 
particle-hole channel effect 

To make this paper self-contained and to introduce the as- 
sumptions as well as the notations, we start by summarizing 
the pairing fluctuation theory [21, 22] without the effect of the 
particle-hole channel, as a foundation for adding the particle- 
hole channel. 

We consider a Fermi gas with a short range s-wave interac- 
tion U (k, k') = U < 0, which exists only between opposite 
spins. Our theory can be effectively represented by a T-matrix 
approximation, shown diagrammatically in Fig. 1 . However, 
we emphasize that this is not a diagrammatic approach, since 
Fig. 1 is simply a representation of the equations derived from 
an equation of motion approach [22, 76-78]. The self energy 
S(i4r) comes from two contributions, associated with the su- 
perfluid condensate and finite momentum pairs, respectively, 
given by T,{K) = T,^{K) + Spg(ii') , where 



Ssc(i^) 



Q 



h^{Q)Go{Q - K) 



(1) 

with Asc being the superfluid order parameter. We use a four- 
vector notation, ii' = (k,ia;(), Q = (q, iri„), Y,k = ^Sz.k' 
etc., and and n„ are odd and even Matsubara frequencies 
for fermions and bosons, respectively. Here Gq{K) = (i — 
$k)"^ and G{K) = [G^^ - T.{K)]~^ are the bare and full 
Green's functions, respectively, ^k = ti^k"^ /2m — jiis the free 
fermion dispersion, measured with respect to the Fermi level. 
In what follows, we wiU set A;b = h = 1. The pseudogap 
T-matrix 



U 



*pg(^) - 1 + UxiQ) 



(2) 



can be regarded as the renormalized pairing interaction with 
pair momentum Q, where 



x{Q) = Y.G{K)Go{Q-K) 



(3) 



K 



is the pair susceptibility. We emphasize that this asymmetric 
form of xiQ) is not an ad hoc choice but rather a natural re- 
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suit of the equation of motion method. The bare Green's func- 
tion Go comes from the inversion of the operator Gq^ which 



appears on the left hand side of the equations of motion. It 
also appears in the particle-hole susceptibility Xph, as will be 
shown below. Albeit not a phi-derivable theory, the equation 
of motion method ensures that this theory is more consistent 
with the Hamiltonian than a phi-derivable theory. 

The gap equation is given by the pairing instability condi- 
tion. 



1 + Ux{0) =0, (T < T,), 



(4) 



referred to as the Thouless criterion, which can also be nat- 
urally interpreted as the Bose condensation condition for the 
pairs, since 1 + Ux{0) oc Upw- In fact, after analytical contin- 
uation iO„ J7 + iO+, one can Taylor expand the (inverse) 
T-matrix as 



ipg(r2,q) 



z- 



(5) 



q^/2M*tothe lead- 



and thus extract the pair dispersion Qq 
ing order, where M* is the effective pair mass. Here Fq is 
the imaginary part of the pair dispersion and can be neglected 
when pairs become (meta)stable [21, 22, 78]. In the superfluid 
phase, tpg{Q) diverges at Q = and a macroscopic occupa- 
tion of the (5 = Cooper pairs, i.e., the condensate, appears. 
This macroscopic occupation, has to be expressed as a singu- 
lar term, tsc(Q) = -(A2^/T)(5(g), (the dashed line in Fig. 1), 
such that ^sc{K) = X^q isc((9)Go((5 - K), written in the 
same form as its pseudogap counterpart, I]pg{K). 

Now we split Y.pg{K) into coherent and incoherent parts: 



Spg(i^) 



E 

Q 



ipg(Q) 



iw; - ^q-k 



E 



^pg(Q) 



ST. 



A2 

pg 



iw; + ^k 



where we have defined the pseudogap Apg via 



A2 
pg 



E^. 

Q 



pg 



ST, (6) 



(7) 



where b{x) is the Bose distribution function. Below T^, the 
divergence of tpg(Q = 0) makes it a good mathematical ap- 
proximation to neglect the incoherent term ST so that 



ioji + ^k 



with 



A^ = Al 



pg ■ 



(8) 



Therefore, the Green's function G{K), the quasiparticle dis- 
persion — y^^^ + A^, and the gap equation, as expanded 
from Eq. (4), follow the same BCS form, except that the total 
gap A now contains both contributions from the order param- 
eter A jc and the pseudogap A^;,. 

For a contact potential, we get rid of the interaction U 
in favor of the scattering length a via m/Ana — 1/C/ + 
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Figure 2. Feynman diagrams for the particle-hole susceptibility Xph 
in the presence of self-energy feedback effect. Panel (b) is identical 
to panel (a), with twisted external legs. The total particle-hole mo- 
mentum P in (a) is equal \.o K -\- K' — Q m (b), with Q being the 
particle-particle pair momentum. 



Ek(l/2ek), where eu 
be written as 



12m,. Then the gap equation can 



1 - 2/(i^J 



1 



(9) 



where ](x) is the Fermi distribution function. In addition, we 
have the particle number constraint, n~2 G{K), i.e.. 



-2E 



^k 



(10) 



where = (1 — Ck/-£']^.)/2 is the BCS coherence factor 

Equations (9), (10), and (7) form a closed set. For given 
interaction l/fcpa, they can be used to solve self consistently 
for Tc as well as A and p at Tc, or for A, Asc, Apg, and p 
as a function of T below Tc. Here fcp is the Fermi wave vec- 
tor. More details about the Taylor expansion of the inverse T 
matrix can be found in Refs. [78, 79]. 



B. Dynamic structure of particle-hole susceptibility Xph(f') 

Before we derive the theory with full particle-hole T-matrix 
tph included, we first study the dynamic structure of the 
particle-hole susceptibility Xph(^')- It is the single rung of the 
particle-hole scattering ladder diagrams, as shown in Fig. 2(a). 
Note that direct interaction exists only between fermions of 
opposite spins. Therefore, the particle and hole must also 
have opposite spins. The total particle-hole four-momentum 
is given by P = {wn, p)- Since we are considering the effect 
on pairing induced by the particle-hole channel, we can twist 
external legs of the diagram, as shown in Fig. 2(b), so that the 
particle-hole contribution can be added to the original pairing 
interaction U directly. It is obvious that the particle-hole mo- 
mentum P in Fig. 2(a) are equal to K + K' — Q in Fig. 2(b), 
where Q is the pair momentum of the particle-particle chan- 
nel. Therefore, we have 



Xph 



K 



E 



/(gk)-/(a-p) ^2 _ i-./(gk)-/(a-p) ^ 

- Ck-p - Wn ^ E^^+ ^k-p + iVn 



(11) 
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Note that again we have a mixing of dressed and undressed 
Green's function in Xvh{P)' like in the expression of x{Q)- 
As mentioned earlier, this mixing has exactly the same origin 
in both cases [80]. For convenience, here we dress the particle 
propagator with self energy and leave the hole propagator un- 
dressed. This is based on the fact that the hole propagator is 
undressed in x{Q) in Sec. II A. (One can equivalently dressed 
the hole while leaving the particle undressed). 

A few remarks are in order. Firstly, in the center-of- 
mass (COM) reference frame of a particle-particle pair, the 
momenta in Fig. 2(b) will be relabeled as ±K + Q/2 and 
±K' + Q/2, so that P = {K + Q/2) - {-K' + Q/2) = 
K + K', independent of the total pair momentum Q. Here 
±K and ±K' are the four momenta of the incoming and 
outgoing fermions in the COM reference frame. Therefore, 
the induced interactions still conform to the Galileo trans- 
formation. Secondly, as in the Nozieres and Schmitt-Rink 
(NSR) theory |3|, one needs a fictitious separable potential 
t^k,k' = U fiifii' in order to have a simple result in the form 
of Eq. (2) for the summation of the particle-particle ladder 
diagrams [17, 21, 22]. The contact potential considered for 
atomic Fermi gases automatically satisfies this requirement. 
However, since the particle-hole pair susceptibility Xph{P) 
only depends on the sum P = K + K',itis obvious that when 
the particle-hole contribution is included, the total effective in- 
teraction ?7eff (k, k') will no longer be separable. Approximate 
treatment is needed, as will be shown later. 

To proceed, we separate the retarded x^ into real and imag- 
inary parts, Xphi'^, P) = XphC^^, P) + iXph(i^- P). after analyti- 



cal continuation, iz^„ 



10+. We have 



XphW, p) 

= E 



/(i?k) - /(Ck-p) ^,2 1 - .f (^k) - /(Ck-p) 
k 



- ?k-p - 



- PJ 2 





+ Ck-p + V 



(12a) 



Xph(^> P) 

= TT^ { - /(i?k - ^)] <S{E^ " Ck-p - ^) 

k 

+ [f{E^ + u)- f{E^)] vl5{E^ + a_p + z.)} . 

(12b) 

From Eq. (12b), we can immediately conclude Xph(O.p) = 
0. In addition, x/i^{v,(^) = if -min(£;k + &) = 
-(-\//"^ + A2 -/i) < V < Q ov V > mayi{E^ - ^k) = 
-\//x^"+A^ -|- fjb. At low T, we also have XphC*^) 0) exponen- 
tially small for \v\ < A if /i > or for \v\ < ^/JjF^TA^ 
otherwise. In all cases, X^i'^^P) = X^(^)P) is isotropic in 
p. In the BCS limit, A ^ 0, ^ |^k|, so that 



x;h(o,P^o)«^/'(a 



"2^ 



(13) 



where /c^ = y/2mjl is the momentum on the Fermi surface. 
For comparison, we also study the dynamic structure of the 



undressed particle-hole pair susceptibility, 

Xp°h(P) = J2 Go{K)Go{K ~P) = Y1 



K 



Ck - ^-p - 

(14) 

This is the particle-hole pair susceptibility studied by GMB 
[61] and others [63-66] in the Uterature. The imaginary part 
is given by 



X, 



o^^^, P) = ^Yl [/(^k) - /(^k - 1^)] <5(a - a-p - ■ 



(15) 

Same as in the Xph(-P) case, we have Xph (^i p) = 0. For finite 
v ^ 0, Xph ('^' P) ^ exponentially as p ^ 0. The real part 
satisfies 



and 



X^Uz/^0,0)=0 



Xp>,P^O)=^/'(a) 



(16) 



2^ 



(17) 



at low T. More generally, for = and finite p, we have 



X%{0,P) = / 
Jo 



kdk m 
2^7 



2k 



2k +p 



(18) 



In the weak coupling limit, Xph(0,P — )• 0) = XpU^'P 0)' 
since Xph reduces to Xph when the gap A vanishes. 

Let's take a look at the case of small but finite p. Equation 
(15) can be rewritten as 



xp°;;>,p) = 7r^[/(a)-/(Ck 

k 



m 



2m 



kp 



m 



P 

2m 



v ] oc v . 



(19) 



The delta function can be satisfied only at fc > \iy\m/p. When 
\v\m/p > fcju, a > for all k so that the magnitude of 
Xph(i^, p) will also turn around and start to decrease expo- 
nentially. The turning points u = ipk^/m show up as two 
peaks in Xph (i^)P)- 

It is easy to show that the hermitian conjugate Xpif * ('^i p) = 
Xpif (— i^, p). Similar relations do not hold for Xph, however, 
due to the mixing of Go and G in the expression of Xph(-f')- 

Shown in Figs. 3(a) and 3(b) are three dimensional (3D) 
plots of the real and imaginary parts of Xpif (i^!P)- In 
Figs. 3(c) through 3(f) we present our calculated results of 
Xj^(j^, p) in the presence of self energy feedback. They are 
calculated at (a-d) and O.lTc (e,f) in the unitary limit, 
1/kpa = 0. Here T^/E^ « 0.256 is the one calculated in 
the pairing fluctuation theory without including the particle- 
hole channel contribution. The even and odd symmetries of 
Xph(i', p) and Xph P) with respect to — )■ — are evident. 
And indeed these symmetries are not present for X^C'^i P)- 
The interesting structure at low frequency and low momentum 
clearly derives from the pseudogap already present at in the 
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Figure 3. (Color online) 3D plots of the real (a,c,e) and imaginary (b,d,f) parts of the particle-hole pair susceptibility Xph with (c-f) and without 
(a,b) the self-energy feedback, calculated at Tc (a-d) and O.lTc (e,f) in the unitary limit. Here Tc is calculated using the pairing fluctuation 
theory without the particle-hole channel contributions. While the undressed Xpif (^iP) has a simple symmetry under v — !> —v, the dynamic 
structure of X^{i',p) is much more complex. In both cases, the real and imaginary parts have very strong dependencies on the total frequency u 
and total momentum p. X^i^^P) shows strong gap effects both at Tc and low T. In units of Bp, the parameters are: Tc = 0.256, /i(Tc) = 0.62, 
A(rc.) = 0.64, ^l{Q.lT,) = 0.59, and A(O.irc) = 0.69. 
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Figure 4. (Color online) Strong momentum dependence of the real 
part of the particle-hole susceptibility at zero frequency v — 
in the unitarity limit, with (black curve) and without (blue dashed 
curve) self-energy feedback, calculated at T = O.STc, where Tc = 
0.256£f. While the undressed Xph(0,p) = ReXph(0,p) shows a 
simple monotonic behavior, the dressed susceptibility Xph(0,p) = 
ReXph(0,p) has a nonmonotonic p dependence, and a substantially 
reduced value at p = 0. This reduction derives from the gap effect 
in the Green's function G{K). Namely, Xph(OiP) seriously over- 
estimated particle-hole fluctuations. 



fully dressed Green's function. In other words, by neglecting 
the feedback effect, the bare Xphi^) niisses this important dy- 
namic structure. The plots of O.IT^ (not shown) is very 
similar to its T = counterpart shown in Figs. 3(a) and 3(b), 
except that the peaks become sharper Comparing Figs. 3(e) 
and 3(f) with 3(c) and 3(d), the gap induced structures become 
much more pronounced. For example, for p — 0, the range of 
1/ in which Xphi'^' 0) = becomes much wider at low T. 

More quantitatively readable two-dimensional plots are 
presented in the Appendix. There we study in detail the effect 
of temperature and interaction strength on the particle-hole 
susceptibility, and how it behaves as a function of frequency 
v for fixed momentum p or as a function of total momentum 
p for fixed frequency 

In Fig. 4, we plot systematically the zero frequency value of 
the real part of the particle-hole pair susceptibiUty as a func- 
tion of total momentum p, with and without the feedback ef- 
fect. The curves are computed at a relatively low T ~ O.STc 
at unitarity. Due to the large excitation gap A = 0.69i?F, at 
p = 0, the value x'phi^' 0) with the feedback is strongly sup- 
pressed from its undressed counterpart, Xp^(0,0). In other 
words, the neglect of the self-energy feedback in xJ^/^(0,0) 
leads to serious over-estimate of the particle-hole channel 
contributions. At the same time, x'-phi^iP) exhibits a more 
complex, nonmonotonic dependence on p than Xpj!|(0,p). In 
both cases, the momentum dependence is strong. 



From Figs. 3 and 4, as well as the Appendix, one can read- 
ily see that the particle-hole susceptibility x^(t^, p) has very 
strong dependencies on both frequency and momentum, as 
well as the temperature and interaction strength. 



(a) it^ = 



O-K O-K' 



(b) M = 



+ 



+ 



-I- 



-I- 



K-P K'-P 

K K 



(c) 



+ 



Q-K Q-K- 



Figure 5. (Color online) Feynman diagrams showing the particle- 
hole effect on fermion pairing, in the presence of self-energy feed- 
back, (a) Particle-particle T matrix ti{Q), with external four mo- 
menta labeled, (b) Particle-hole T matrix tph(P), with P = K + 
K' — Q being the total particle-hole 4-momentum. (c) An effec- 
tive, composite particle-particle T-matrix, t2{Q), with the contribu- 
tion from the particle-hole channel included. Here different shadings 
represent different T matrices. 



C. Induced interaction - beyond the lowest order 

Figure 3(b), in the absence of the self-energy feedback ef- 
fect, is in fact the lowest order induced interaction, considered 
in GMB and most others in the literature: 

uUp) - -uXuiP)- (20) 

Diagrams of the same order but between fermions of the same 
spin vanish. 

Let us first re-plot the particle-particle scattering T-matrix 
in Fig. 5(a); this is just tpg in Fig. 1 but with external legs, 
which we here refer to as ti {Q). We have 



1 



xiQ) 



(21) 



Now we consider the contribution of an infinite particle-hole 
ladder series, as shown in Fig. 5(b). Such contribution should 
be added to the bare interaction U. The summation of the se- 
ries of particle-hole ladder diagrams gives rise to the T-matrix 
in the particle hole channel. 



u 



1 



(22) 



1 + UxphiP) + Xph(P) 

At Q = 0, this gives the overall effective pairing interaction, 

U,niK, K') - tp^iK+K') = U+U^ ^ 



1 + Uxph{K + K') ' 
(23) 

where K and K' are the incoming and outgoing 4-momenta 
of the scattering particles in the COM reference frame. The 
induced interaction is thus given by 



C/i„d(T) = tph(T) - U 



1 + C/Xph(P) ' 



(24) 
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with P = K + K' . Upon Taylor expanding the denomina- 
tor by treating ?7xph as a small parameter, one readily notice 
that the leading term is just the numerator This would be the 
counterpart lowest induced interaction in our theory, except 
that we always consider the self energy feedback effect. 

It is evident that the T matrices in the particle-particle chan- 
nel and the particle-hole channel share the same lowest order 
term, U. Both T matrices can be regarded as a renormalized 
interaction, but in different channels. What we need is to re- 
place the bare U in one of the two T matrices with the other 
T matrix. The results are identical, which we call t2. Shown 
in Fig. 5(c) is the regular particle-particle channel T matrix 
ti{Q) with U replaced by the particle-hole channel T matrix 
fph(P) (with twisted external legs), where P = K + K' — Q. 
In other words, we replace 'witht~^{P) = U~^+Xph{P) 
in Eq. (21), and formally obtain 



t2{Q) = 



1 



U-^+xAK + K'-Q)+x{Q) 



(25) 



Unfortunately, since Ueff{K, K') is not a separable potential, 
one cannot obtain a simple summation in the form of Eq. (25). 
This can also be seen from the extra dependence on K and 
K' on the right hand side of the equation. Certain averaging 
process has to be done to arrive at such a simple summation, 
as will be shown below. 



D. Gap equation from self-consistency condition in mean-field 
treatment 

The dependence of Uef[{P) on external momenta via P = 
K + K' — Q presents a complication in the gap equation. 
This can be seen through the self consistency condition in the 
mean field treatment, even though we do not use mean field 
treatment in our calculations. Writing the interaction Vk,k' = 
Uee{K + K') for (5 = 0, i.e., zero total pair 4-momentum, we 
have 



K' 



Vk,K'{Ck'C_k') 



u 



^l + UxAK + K'){iuJvY-El, 



, (26) 



where we have used the mean-field result {cj(,c_j(,) = 
G{K')Go {—K')Ak' . Equivalently, this can be written as 



K' 



u 



A. 



+ [/xph(iwi +iw;/,k + k') {iuj'f 



(27) 

Note that, due to the dynamic character of Xph(-fi' + K'), both 
the gap Ax and the quasiparticle dispersion Ek acquire a dy- 
namical frequency dependence. The gap also develops a mo- 
mentum dependence, which is originally absent for a contact 
potential. 

We can express Uqs{P) in terms of its retarded analytical 
continuation, as follows: 



dv -2Imt/Jf(i/,p) 



(28) 



where the second term is just the induced interaction, 

Xph(^^'P) 



Im?7e^(z/,p) = 



{u-' + x!,^? + {x!^?- 



(29) 



Then we have 



ILJi + IW;/ - 

(30) 

The particle-hole channel effect is contained in the 2nd term, 
without which this would be just the gap equation without 
the particle-hole channel, and admit a constant gap solution. 
Without further approximation, the complex dynamic struc- 
ture of Xph{P) will inevitably render it very difficult to solve 
the gap equation. 



E. Pairing instability condition in tlie presence of particle-hole 
channel effect 



In order to obtain a simple form as Eq. (25), we have to 
average out the dependence of Uss{K, K') on K and K' . In- 
deed, an average of Xph(i^, p) has been performed in the Uter- 
ature on (and only on) the Fermi surface |61|. For the fre- 
quency part, here we follow the literature and take ivn = 
ioji + iuif = 0. From Fig. 3, one can see that this is where 
the imaginary part Xphi'^^P) — ^ P ^h^^ the effec- 

tive interaction UsifiK, K') is purely real. For the momentum 
part, we choose on-shell, elastic scattering, i.e., k = k', and 
then average over scattering angles: 



p = |k + k'l = fcA/2(l + C0S6') , 



(31) 



where 9 is the angle between k and k'. It is the off-shell scat- 
tering processes which lead to imaginary part and nontrivial 
frequency dependence in Xpui^^P) order parameter. 

Further setting k = k^ and averaging only on the Fermi sur- 
face is the averaging process used in all papers we can find 
about induced interactions in the literature. We refer to this 
as level 1 averaging. In this paper, we also perform a level 
2 average, over a range of k such that the quasiparticle en- 
ergy E]^ € [min(i!^i^), miii(£'^) + A]. Here mm{E^J = A 
if /X > 0, or min{E^^) = v'm^ + A^ if fi < 0. The ba- 
sic idea is that according to the density of states of a typi- 
cal s-wave superconductor, the states within the energy range 
Ej^ e [A, 2A] are most strongly modified by pairing. It 
should be pointed out that in the BEC regime, this range can 
become very large. 

Upon averaging of either level 1 or level 2, we drop out 
the complicated dynamical structure of Xphii^jP) and replace 
it by a constant (xph)- For the purpose of comparison, we 
shall also perform the averaging on the undressed particle- 
hole susceptibility Xphi'^' P) mostly show the result at 
level 1. 

Shown in Fig. 6 are the angular averages of the particle- 
hole susceptibility at = as a function of momentum k 
under the above on-shell condition, k = k'. Here we only 
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\lk^ = 




<x^^>,Trr, = o.i 

<XpV,7-/7-, = 0.1 



k/k. 



Figure 6. (Color online) Angular average of the on-shell particle- 
hole susceptibility, (xph(0,p = |k + k'|)) at = as a function 
of momentum k/kf, under the condition k — k' , calculated at uni- 
tarity for different temperatures T — O.lTc (black solid curve) and 
T — (green dot-dashed curve), in units of k^/Ej. Also plotted 
is its undressed counterpart, (Xph(0,P = |k + k'|)), which shows a 
serious over-estimate due to the neglect of the self-energy feedback. 
Here Tc = 0.256i5F and associated gap and /i values are calculated 
without the particle-hole channel effect. The open circles on each 
curve denote level 1 average, i.e., k — k^. The vertical axis readings 
of the horizontal short bars indicate the corresponding values of level 
2 average. The thick section of each curve indicates the range of k 
used for level 2 averaging. Clearly, there are strong temperature and 
k dependencies in both (Xph(0,p)) and {Xph(0,p)). The (absolute) 
values of Level 2 average are substantially smaller than their level 1 
counterpart. 



show the unitary case at two different temperatures, T — 
and low T = Q.IT^ ^ T^. For the purpose of comparison, 
we plot the result for both the dressed and undressed particle- 
hole susceptibility. The curves show strong momentum de- 
pendencies. For (Xph('^'P))' i'^ monotonically increasing, 
whereas for (xph(0,p)), it exhibits nonmonotonic k depen- 
dence at low T. Both dressed and undressed particle-hole sus- 
ceptibilities have a temperature dependence, and this depen- 
dence is much stronger for the former. This can be attributed 
mainly to the temperature dependence of A(T) in (xph(0,p)), 
while (Xpjj(0,p)) depends on T only via 

The open circles on each curve represent the level 1 aver- 
age, i.e., the values at fc = A;^. At the same time, the ver- 
tical axis readings of the short horizontal bars correspond to 
the level 2 average, while the thick segments of each curve 
represents the range of k used for level 2 averaging. Figure 
6 shows that the (absolute) values of the level 2 average are 
significantly smaller than their level 1 counterpart. The level 
1 average {Xph{0,p)) is essentially temperature independent 
(see the red and blue circles). In addition, it is evident that the 
neglect of self-energy feedback has caused {x^h i^-iP)) -^e'"'- 
ously over-estimate the contribution of particle-hole channel. 

Similar plot for l/fcpa = 0.5 (Fig. 15 in the Appendix) 
exhibits a much stronger T dependence. In that case, fi is 
very close to albeit still positive. As a consequence, the 
particle-hole susceptibility is much smaller than that shown in 



Fig. 6. 

Now with this frequency and momentum independent 
Xph(j^,p) ~ (Xph), we can easily carry out the simple geo- 
metric summation for 12'. 



1 



u-^ + (xph) + xiQ) 



(32) 



Therefore, the Thouless criterion for pairing instability leads 
to the gap equation: 



+ (Xph) + X(0) = 0, 



namely. 



/ m 
V47ra 



(xph)) - E 



1 - 2/(^1,) 1 
2E^ 2ek 



(33) 



(34) 



As will be shown later, (xph) is always negative. Therefore, 
the particle-hole channel effectively reduces the strength of 
the pairing interaction. 

In the weak coupling limit (l/fcpa = — oo), A ^ 0, T < 
Tc ^ Tp, then (xph) and (Xph) become equal, for either level 
of averaging. We have 



(Xph) = 



J J' Jo i^7^(^''^^" 



2k -p 



2k +p 



N{0) 





2k 


-P 


f dx f In 




Li Jo 2p 


2k 


+ P 



1±|^A.(0) = 0.02015# 
3 Ep 



p=kr^2 (1+x) 



(35) 



where k ~ k/kp, p — p/kp ~ \^2{1 + x), x — cos 9, and 
A^(0) — mkp/2TT'^h^ is the density of state at the Fermi level. 
Here we have approximated the Fermi function with its T = 
counterpart, with a step function jump at the Fermi level. 

In the weak interaction limit, the BCS result for is 
T^^^Ep = (8/7r)e^-2gi/Ar(o)c/^ where 7 w 0.5772157 is 
the Euler's constant. Equation (33) implies a replacement of 
1/f/by 1/U + (xph). In this way, the new transition temper- 
ature Tc is given by 



.{Xph>/A^(0) 



(4e) 



-1/3 



J'BCS 

and the same relation holds for zero T gap, 
A 



0.45, 



= (4e) 



-1/3 



(36) 



(37) 



This result is in quantitative agreement (to the leading order) 
with that of GMB [61] and others [63] in the literature. Note 
that in our work, as well as in that of Yin and coworkers [66], 
the average particle-hole susceptibility (xpti) is added to 1/C/ 
or m/47ra. In other works [63-65], only the lowest order 
particle-hole diagram is considered so that their induced in- 
= —W^iXph) is added to U. Therefore, these 



teraction U^^^ 



works have to rely on the assumption N{0)U <^ 1 and the va- 
lidity of the BCS mean-field result in order to obtain the result 
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of Eq. (36). Away from the weak interaction regime, a full 
summation of the particle-hole T matrix becomes necessary. 

While the results from all different treatment seem to agree 
quantitatively in the weak coupling limit, we expect to see 
significant departures as the pairing interaction strength in- 
creases, especially in the unitary regime. 

With the overall effective interaction C/eff, the self energy, 
as obtained from T.{K) = J2q h{Q)Go{Q - K), will follow 
the same form as Eq. (8) although the gap values will be dif- 
ferent. Therefore, the fermion number equation will also take 
the same form as Eq. (10). Furthermore, the pseudogap equa- 
tion, given by Ap^ = — J^q t2{Q), will also take the same 
form as Eq. (7). 

Equations (10), (7), and (34) now form a new closed set, 
and will be solved to investigate the effect of the particle-hole 
channel. 

Note that in a very dilute Fermi gas shifting m/ATra by 
(Xph) has no significant influence in experimental measure- 
ment of the s-wave scattering length a, because (xph) has di- 
mension [fcF]'^/[£^F] — [kp] and thus vanishes as fcp — > in 
the zero density limit. However, a finite kp will indeed shift 
the resonance location except at very high T where /i < 0. In 
Ref. [56], from which the scattering lengths are often quoted 
for ^Li, the actual density is comparable or even higher than 
that in most typical Fermi gas experiments. Therefore, the 
particle-hole channel may play an important role. 

Here we propose that this particle-hole channel effect may 
be verified experimentally by precision measurement of the 
magnetic field B at the exact Feshbach resonance point as a 
function of density or kp at low T. The zero density field So 
can be obtained by extrapolation. Then one should have the 
field detuning SB — B — Bq oc kp. Because different the- 
ories predict a very different value of (xph) at unitarity, the 
measured field detuning can thus be used to quantify (xph) 
and test these theories. In principle, one may experimentally 
measure (xph) through the entire BCS-BEC crossover. For a 
Fermi gas in a trap, the trap inhomogeneity leads to a distri- 
bution of kp. Instead of a uniform shift, this inhomogeneity 
will spread out the unitary point at zero density into a narrow 
band at finite density. The band width and mean shift are both 
expected to be proportional to kp. Such effect deserves further 
investigation. 



III. NUMERICAL RESULTS AND DISCUSSIONS 
A. Effect of particle-hole channel on BCS-BEC crossover 

In this section, we will investigate the effect of the particle- 
hole channel on the BCS-BEC crossover behavior, in terms of 
zero temperature gap A(0), and their ratio. 

First, in Fig. 7, we show the effect on the zero T gap by 
comparing the calculated result with and without the particle- 
hole channel contributions. Shown respectively in panel (a) 
and (b) are plots of the zero T gap A and the corresponding 
particle-hole susceptibility (with a minus sign) as a function 
of 1/kpa. The black solid line in Fig. 7(a) is the result with- 
out the particle-hole channel effect, whereas the other curves 
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Figure 7. (Color online) Effect of the particle-hole channel contri- 
butions on the zero temperature gap in BCS-BEC crossover. In (a), 
the black solid curve is the gap without the particle-hole effect. The 
rest curves are calculated with the particle-hole channel effect but 
at different levels, i.e., using undressed particle-hole susceptibility 
(Xph) with level 1 averaging (red dotted line), dressed particle-hole 
susceptibility (xph) with level 1 (green dot-dashed curve) and level 2 
(blue dashed line) averaging, respectively. The corresponding values 
of the average particle-hole susceptibility with a minus sign are plot- 
ted in (b), in units of kp/Ep. The particle-hole channel effect can be 
essentially neglected beyond l/fcpa > 1.5. 



are calculated with the effect at different levels of approxi- 
mation. The (red) dotted curve are calculated using the un- 
dressed susceptibility (Xph) ^'^ average level 1. The (green) 
dot-dashed and (blue) dashed curves are calculated using the 
dressed particle-hole susceptibility (Xph) with level 1 (green 
dot-dashed curve) and level 2 (blue dashed line) averaging, re- 
spectively. The level 2 result shows a slightly weaker particle- 
hole channel effect, as can be expected from Fig. 6. 

One feature that is easy to spot is that the undressed 
particle-hole susceptibility (Xph) has a very abrupt shut-off 
where the chemical potential fi changes sign. As a result, 
the corresponding (red dotted) curve of the gap also merges 
abruptly with the (black solid) gap curve calculated without 
particle-hole channel effect. This is not unexpected as one 
can see from Eq. (18) that (x°h> = at T = for ^ < 0. 
Furthermore, Eq. (17) implies that (Xph) approaches zero at 
/i = abruptly with a finite slope as k^ does. In contrast, 
with the self-energy feedback included, either level 1 (green 
dot-dashed curves) or level 2 (blue dashed curves) average of 
(Xph) approaches smoothly as the BEC regime is reached. 
Consequently, in Fig. 7(a), the (green) dot-dashed and (blue) 
dashed curves approach the (black) solid curve very gradually. 
It is also worth pointing out that the difference between level 1 
and level 2 average of (xph) is less dramatic than that between 
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(Xph) and the undressed (Xph)- Indeed, the (green) dot-dashed 
and (blue) dashed curves are very close to each other. 

The abrupt shut-off of (Xph) at /i = is determined by 
the step function characteristic of the Fermi function at T = 
0. At finite T, this shut-off will become smoother with an 
exponential tail on the BEC side. 

In the unitary regime, especially for 1/kpa E [—0.5, +0.5], 
the particle-hole susceptibility is strongly over-estimated by 
the undressed (Xph) comparison with the dressed (xph)- In 
this regime, both the gap and the underlying Fermi surface (as 
defined by the chemical potential) are large, so that neglect- 
ing the self-energy feedback leads to a strong over-estimate of 
(Xph)' because the large gap serves to suppress particle-hole 
fluctuations. 

From Fig. 7, we conclude that the particle-hole effect di- 
minishes quickly as the Fermi gas is tuned into the BEC 
regime with increasing pairing interaction strength. Beyond 
l/kfG > 1.5, the effect can essentially be neglected. For the 
level 1 average of the undressed particle-hole susceptibility, 
(Xph)' been done in the literature, this effect disap- 

pears immediately once the BEC regime (defined by /i < 0) 
is reached, as far as the zero T gap is concerned. 

As a consistency check, we notice that in the BCS limit, 
the average particle-hole susceptibility in all cases in Fig. 7(b) 
approaches the same asymptote, which is given by Eq. (35). 
This confirms our previous analytical analysis. 

Next, we show in Fig. 8 the effect of the particle-hole chan- 
nel on the behavior of as well as the pseudogap at T^. Fig- 
ure 8(c) can be compared with Fig. 7(b). The curves for levels 
1 and 2 average of (xph) in Fig- 8(c) are very similar to those 
in Fig. 7(b), with the values at 1/kpa = slightly smaller. 
On the other hand, the curve for (Xph) ^ smooth thermal 
exponential tail in the BEC regime in Fig. 8(c). Thus, the 
pseudogap A(Tc) calculated using (Xph) now merges back to 
the (black) solid curve smoothly. 

Similar to the zero T gap case in Fig. 7, the difference in 
the effect on Tc and A(Tc) between level 1 and level 2 averag- 
ing mainly resides in the unitary regime, and is less dramatic 
than that between undressed and dressed particle-hole suscep- 
tibility. Again, the undressed particle-hole susceptibility gives 
rise to an overestimate of the particle-hole channel effect. 

In all cases, the particle-hole susceptibility becomes negli- 
gible in the BEC regime. The effect of the particle-hole chan- 
nel shifts the Tc and A(rc) curve towards larger l/fcpa, al- 
though the amount of shift clearly depends on the value of 
l/fcpa. 

Now we study the effect of the particle-hole channel on the 
ratio 2A{0)/Tc. It suffices to consider the mean-field ratio, 
2A(0)/T;'^^ since 2A(0)/r, obviously wifl deviate from the 
weak coupling BCS result when pairing fluctuations are in- 
cluded in the crossover and BEC regimes. From Fig. 6, we 
see a strong T dependence of the particle-hole susceptibility. 
Therefore, the effect on T^^ and on zero T gap A(0) are dif- 
ferent, as can be seen roughly from Figs. 7 and 8. 

In Fig. 9, we plot this mean-field ratio as a function of 
l/fcpfl with (black solid curve) and without (blue dashed 
curve) the particle-hole channel effect. Here the particle-hole 
susceptibility (xph) is calculated with level 2 averaging. In the 
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Figure 8. (Color online) Effect of the particle-hole channel contribu- 
tions on Tc and the pseudogap A at Tc in BCS-BEC crossover. In (a) 
and (b), the black solid curves are calculated without the particle-hole 
effect. The rest curves are calculated with the particle-hole channel 
effect but at different levels, using undressed particle-hole suscepti- 
bility (Xph) with level 1 averaging (red dotted line), dressed particle- 
hole susceptibility (xph) with level 1 (green dot-dashed curve) and 
level 2 (blue dashed line) averaging, respectively. The corresponding 
values of the average particle-hole susceptibility with a minus sign 
are plotted in (c), in units of kf /Ep. The particle-hole channel effect 
can be essentially neglected beyond l/kpa > 1.5. 



l/kpa — oo limit, the ratio is unaffected by the particle- 
hole channel. As l/kpa increases, the contribution of the 
particle-hole channel causes this ratio to increase gradually. 
At l/kpa = —4, which is still a very weak pairing case, the 
ratio is already slightly larger. The effect is most dramatic in 
the unitary regime, since further into the BEC regime, (xph) 
will vanish gradually. It is worth noting that even without the 
particle-hole channel, the ratio 2A(0) /T^^ starts to decrease 
from its weak coupling limit, 2Tre^'' w 3.53. 

Finally, we estimate the shift in Feshbach resonance posi- 
tions. From Figs. 7 and 8, we find that Xph does not neces- 
sarily diminish as T increases except at very high T (where 
/i becomes negative, so that |xph| will decrease exponentially 
with T.) In fact, this can be understood because A(r) de- 
creases with T so that |xph| increases. We take (xph) = 
-O.Olkl/Ep = -0.01{2mkp/h'^). According to Eq. (34), 
the shift in 1/a is S{l/a) = -47r/i^(xph)/2TO = COSTrfep- 
In other words, the dimensionless shift S{l/kpa) — 0.25, 
which is independent of density and is no longer negligi- 
ble. This is in good agreement with the actual shift 0.32 of 
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Figure 9. (Color online) Effect of the particle-hole channel contri- 
butions on the ratio 2/^{0)/T^^ in BCS-BEC crossover. Shown is 
the mean-field ratio calculated with (black solid curve) and without 
(blue dashed curve) the particle-hole channel contributions. Here the 
particle-hole susceptibility (xph) is calculated with level 2 averaging. 



the peak location of the Tc curve in Fig. 8(a). For a typ- 
ical Tp = l/^K in ^Li, using the approximate expression 
a = abg[l — W/{B — Bq)], we obtain the shift in resonance 
position 5Bo -0.08M^(fcFabg) = 7.8G. Here for the lowest 
two hyperfine states, the resonance position Bq = 834. 15G, 
the resonance width W — 300G, and the background scat- 
tering length flbg = — 1405ao^ with cq = 0.528A. Clearly, 
the shift (5i?o is not small. In reality, one needs to solve self- 
consistently the equation m/ (47ra) + (Xph) — 0, and take care 
of the trap inhomogeneity. These will likely make the actual 
average shift smaller. 

The susceptibility Xph calculated with and without the self 
energy feedback differs by roughly a factor of 2 at unitarity. 
This can be used to test different theories, as mentioned ear- 
lier. 

A question arises naturally as to whether the particle-hole 
channel effect has already been included in the experimentally 
measured scattering length a, since, after all, the measure- 
ments of a such as those in Ref. [56] were carried out at den- 
sities comparable to typical Fermi gas experiments. This also 
depends on whether the temperature was high enough during 
the measurements. 



B. Critical temperature Tc at unitarity 

Finally, we compare our result on the critical superfluid 
transition temperature T^/E-p for a 3D homogeneous Fermi 
gas at unitarity with those reported in the literature. From 
Fig. 8, we read T^/Ep = 0.217 using level 2 average of (xph)- 
And the maximum Tc w 0.257 now occurs at 1/kpa w 0.32, 
on the BEC side. The level 1 average of (Xph) yields a 
slightly lower value, T^/Ep — 0.209. However, we empha- 
size that the level 2 average of (xph) is more reasonable. Note 
that as in the theory without particle-hole channel effect, we 



have dropped out the incoherent part of the self-energy from 
particle-particle scattering. Inclusion of the incoherent part is 
necessary in order to obtain the correct value of the /3 factor. 

Hu and coauthors [81, 82] have been claiming to be able to 
obtain the coiTect value of the /3 factor, using an NSR-based 
approach, without including the particle-hole channel. Obvi- 
ously, their claim will breakdown when the particle-hole chan- 
nel is included. 

The value of Tc for a homogeneous Fermi gas at unitarity 
has been under intensive study over the past few years. With- 
out the particle-channel effect, using the pairing fluctuation 
theory, we previously reported a value of Tc/Ep = 0.255 for 
a short range Gaussian potential (using a two-channel model 
with a cutoff momentum ko/kp = 80) [29] or 0.256 for an 
exact contact potential. In all cases, the maximum Tc oc- 
curs very close to but slightly on the BEC side of the uni- 
tarity. The original NSR theory [3] predicted TjEp = 0.222. 
Haussmann et al. [83] found 0.16, using a conserving ap- 
proximation which involves only dressed Green's functions. 
However, this theory did not contain a pseudogap at Tc and 
exhibited unphysical non-monotonic first order like behavior 
in the temperature dependence of entropy S{T). Floerchinger 
et al. [84] found a high value of 0.264 even after includ- 
ing the particle-hole channel fluctuations. By including the 
induced interaction in an NSR-based treatment without self- 
energy feedback, Yin and coworkers [66] reported a value of 

0. 178. However, as shown above, the neglect of the pseu- 
dogap self energy feedback in Ref. [66] has caused a serious 
over-estimate of the contributions of the induced interaction. 
In addition, their work suffers all defects of the NSR theory, 

1. e., lack of self-consistency and neglect of the pseudogap ef- 
fect. 

Using the renormalization group method, Stoof and 
coworkers [85] obtained T^/Ep = 0.13, lower than most other 
calculations. Troyer and coworkers [86] reported 0.152 using 
quantum Monte Carlo (QMC) simulations for lattice fermions 
at finite densities and then extrapolated to zero density limit. 
Using (and improved upon) the method of Ref. [86], Goulko 
and Wingate [87] found a higher value, 0. 17 1 . An even higher 
value would have been obtained had they used an quadratic fit 
to (and extrapolation of) their intermediate density data. Also 
using QMC, Akkineni et al. [88] recently reported a value of 
0.245. 

Experimentally, only Tc in a trap has been measured. The 
Duke group [71], in collaboration with Chen et al, found 
Tc/Ep = 0.27 in a trap through a thermodynamic measure- 
ment in a unitary ^Li gas. Later, the Duke group [89, 90] ob- 
tained 0.29 and 0.21 by fitting entropy and specific heat data 
with different formulas. At unitarity, it has been known that 
Tc/ Ep in the trap is only slightly higher than its homogeneous 
counterpart. Therefore, these experimental values probably 
imply that the homogeneous value of T^/Ep is about 0.25 ^ 
0.19. Our result, T^/ Ep = 0.217, is in reasonable agreement 
with these experiments. Recently, Ku et al. [91] reported 
Tc/Ep w 0.167 for a homogeneous Fermi gas by identifying 
the lambda-like transition temperature. Interestingly, if we 
take an constant 5Yj = — 0.3£'f (half of the energy of a single 
spin down atom in a spin up Fermi sea [92]) as the incoherent 
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Figure 10. (Color online) Higher order T-matrix, ts, obtained by 
repeating the T-matrix t2- 

part of the self energy in Eq. (6), T^/E^ will be suppressed 
down to 0.174, close to this recent experimental value. Full 
numerical inclusion of (5S will be done in a future study. 



C. Higher order corrections 

In addition to non-ladder diagrams, which we have cho- 
sen not to consider, there seem to be a series of higher or- 
der corrections. For example, one can imagine repeating the 
T-matrix t2 in the way shown in Fig. IIIB, and obtaining a 
higher order T-matrix tj,. Such can then be repeated to 
obtain a higher order T-matrix t4, and so on. While one may 
argue these higher order T-matrices are indeed of higher order 
in bare interaction U, our experience with t2 seems to imply 
that detailed study needs to be carried out before we jump to 
a conclusion on this. Indeed, even the lowest order so-called 
induced interaction [7;° ^ is one order higher in U than U itself. 



regime. 

Without including the incoherent part of the self energy 
from particle-particle scattering, our present result on the crit- 
ical temperature at unitarity yields T^/Ep w 0.217, substan- 
tially lower than that obtained without the particle-hole effect. 
This value agrees reasonably well with some existing experi- 
mental measurement. 

We have also made a falsifiable proposal that the particle- 
hole contribution can be measured by locating the Feshbach 
resonance positions as a function of fcp and that this can be 
used to test different theories. 

To study more accurately the quantitative consequences 
of the dynamic structure of the particle-hole susceptibility, 
full-fledged numerical calculations are needed, without taking 
simple angular average and setting frequency i/ = 0. Further 
investigation is called for in order to determine whether higher 
order T-matrices will make a significant difference or not. 
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IV. CONCLUSIONS 

In summary, we have studied the effects of the particle-hole 
channel on BCS-BEC crossover and compared with lower 
level approximations. We include the self-energy feedback in 
the particle-hole susceptibility Xph, which leads to substantial 
differences than the result without self-energy feedback. 

We have investigated the dynamic structure of Xph, and 
have discovered very strong temperature, momentum and fre- 
quency dependencies. Angular (as well as radial) average in 
the momentum space of the particle-hole susceptibility has 
been done in order to keep the equations manageable. We 
have performed the average at two different levels and also 
compared with the result calculated without including the self- 
energy feedback. We conclude that the level 2 averaging, i.e., 
both over angles and a range of momentum, is more reason- 
able. Computations of the particle-hole susceptibility with- 
out the self-energy feedback leads to an overestimate of the 
particle-hole channel effect. 

In the weak coupling BCS limit, our result agrees, to the 
leading order, with that of GMB and others in the litera- 
ture. Away from the weak coupling limit, A(0) and Tc 
are suppressed differently. We have also studied the ra- 
tio 2A(0)/T^'^^ at the mean-field level and found that it is 
modified by the particle-hole fluctuations. The particle-hole 
channel effects diminish quickly once the system enters BEC 



Appendix A: Dynamic structure of the particle-hole 
susceptibility 

In this Appendix, we will present a series of two- 
dimensional plots, in order to make the 3D data shown in 
Fig. 3 quantitatively easier to read. 

First, we study the impact of temperature. We present in 
Fig. 11 the real and imaginary parts of Xphi^^O) for differ- 
ent T from low to high in the unitary limit, l/fcpa = 0. To 
single out the temperature effect, here we take for all tem- 
perature A ~ 0.686 and /i = 0.59, which are their val- 
ues calculated at T = using the pairing fluctuation the- 
ory without the particle-hole channel effect. Evidently, be- 
yond the point v — ^/^^ + A^ + fi ~ 1.49, the imagi- 
nary part Xphi'^^^) vanishes identically. Around v — 0, the 
lower bound of the range of where Xphi'^^^) essentially 
vanishes changes from u = —A —0.686 at low T to 
ly ^ -(Vm^ + A2 - m) = 0.315 at high T. Meanwhile, 
its upper bound decreases continuously with T from = A at 
low T to 1/ = at very high T. This numerical result agrees 
with our previous analysis. A comparison with the real part re- 
veals that the peaks in Xph(^i 0) correspond to the sharp rises 
in the plot of Xphi'^' 0) ^^^^ these lower and upper bounds. 
This can also be seen from the Kramers-Kronig relation be- 
tween the real and imaginary parts of 0). 
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Figure 11. (Color online) The particle-hole susceptibility Xph('^iO) 
at zero momentum p in the presence of self-energy feedback for dif- 
ferent T (as labeled) at unitarity. To single out the temperature effect, 
we fix A = 0.686 and /i = 0.59 at their values at T = 0, calculated 
using the pairing fluctuation theory without the particle-hole channel 
effect. 



In comparison, we have also studied the temperature evolu- 
tion of the undressed Xpjf (j^, p). Shown in Fig. 12 is the result 
for p = 0.1. From Fig. 3, it is easy to see that one cannot plot 
the result forp ~ 0. This can also be seen from Eqs. (16) and 
(17). For finite p, say p = 0.1, the peaks at low T in both real 
and imaginary parts become more smeared out as T increases. 
Near = 0, we see that Xph (^'P = ^■^) proportional to 
V, in agreement with our previous analysis. Xpif ('^iP) shows 

good symmetry about v: x^hi^^^P) = Xpif'*('^iP)- There is 
no gap effects, of course. 

Shown in Fig. 13 is the evolution of the particle-hole sus- 
ceptibility x^(i^, 0) at total momentum p = in the presence 
of feedback effect with increasing pairing strength. These 
curves are calculated at low temperature T — O.lTc. Here 
for each interaction strength, the parameters A, /i and Tc 
are calculated using the pairing fluctuation theory without 
the particle-hole channel effect. Around v ~ {), the range 
within which the imaginary part vanishes is given by < A 
for the /i > cases (l/fcpa = —1 through 0.5). For 
l/Zcpa = 1, IJ-/Ef = —0.8, the lower bound is given by 
— {y^Jjp + A? — /i) = —2.35 and its upper bound extends 
to oo since A = 1.33 > ^/i^ + A^ + /i = 0.75. It is obvi- 
ous that this range becomes wider and wider with increasing 
pairing strength from BCS to BEC. 

From Figs. 3(a) and 3(b), one readily notice that for the 
undressed Xph'{i',p), it is not appropriate to plot Xphi'^^P — 
0) as a function of i^. Instead, one has to plot at a finite p, say, 
p = 0.1 fcp, in order to study its temperature evolution. Our 
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Figure 12. (Color online) The undressed particle-hole susceptibility 
= 0.1) for different T (as la 
0.1. As in Fig. 11, we fix : 



Xph*('^i P = 0.1) for different T (as labeled) at unitarity and momen^ 



tum ) 



0.59. 



result (not shown) demonstrates that the real part presents a 
double peak structure, with the peaks becoming increasingly 
broader as T increases. At low T, the location of the peaks 
are roughly given by = ±.pk^lm « ±0.15 for /i — 0.59 at 
unitarity. This relation also shows how the ^^{v^p) curves 
evolve with total momentum p. 

Next, we investigate how the particle-hole susceptibility 
X,^(i^, p) evolves with total momentum p in the presence of 
feedback effect. Shown in Fig. 14 are the curves of the real 
and imaginary parts for increasing p for a unitary Fermi gas, 
calculated at Tc. Just as in Fig. 11, the p = curve shows 
a clear gap in the neighborhood of = in the imaginary 
part, x!^\Syi'P)- As p increases, this gap gradually disappears, 
and the upper bound in v beyond which Xph('^iP) vanishes in- 
creases towards infinity. At the same time, the peaks in the 
real part becomes broader and smeared out. From Fig. 14(a), 
we can see that at = 0, the real part slowly increases with p. 

The zero frequency value Xph(0,p) is plotted in Fig. 4 in 
the text as a function of p. 

Finally, we show in Fig. 15 the angular average of the on- 
shell particle-hole susceptibility, (xph(0,p — |k + k'|)) at 
I' = as a function of momentum fc/fcp, under the condi- 
tion k = k', calculated at l/fcpa = 0.5. The chemical poten- 
tial is nearly zero, close to the boundary separating fermionic 
and bosonic regimes. In comparison with the unitary case 
shown in Fig. 6, we conclude that both dressed and undressed 
particle-hole susceptibility exhibit stronger temperature and 
fc dependence. Here the small chemical potential determines 
that the susceptibility is also much smaller. It is worth men- 
tioning that the level 1 average of the undressed particle-hole 
susceptibility actually shows a much stronger temperature de- 
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Figure 13. (Color online) The real and imaginary parts of particle- 
hole susceptibility Xja{^j 0) the presence of self-energy feedback 
for various values of l/fcpa from BCS to BEC. The curves are cal- 
culated at O.lTc. For each case, the parameters A, fi and Tc are cal- 
culated using the pairing fluctuation theory without the particle-hole 
channel effect. 
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Figure 15. (Color online) Angular average of the on-shell particle- 
hole susceptibility, (xph(0,p = |k -1- k'|)) at v — as a func- 
tion of momentum fc/fcp, under the condition k — k' , calculated at 
l/kpa = 0.5. The conventions and legends are the same as in Fig. 6. 
Here Tc = 0.226£'f and the chemical potential /i/Bp = 0.077 and 
0.074 at Tc and O.lTc, respectively. Clearly, there are even stronger 
temperature and k dependencies in both (Xph(0,p)) and {Xph(0,p)) 
than the unitary case shown in Fig. 6. The (absolute) values of Level 
2 average are substantially smaller than their level 1 counterpart. 



16 



[69] T. Timusk and B. Statt, Rep. Prog. Phys. 62, 61 (1999). 
[70] R. Damascelli, Z. Hussain, and Z.-X. Shen, Rev. Mod. Phys. 
75, 473 (2003). 

[71] J. Kinast, A. Turlapov, J. E. Thomas, Q. J. Chen, J. Stajic, and 
K. Levin, Science 307, 1296 (2005), published online 27 Jan- 
uary 2005; doi: 10. 1126/science. 1109220. 

[72] Y. He, Q. J. Chen, and K. Levin, Phys. Rev. A 72, 011602(R) 

(2005). 

[73] Q. J. Chen and K. Levin, Phys. Rev. Lett. 102, 190402 (2009). 
[74] J. T. Stewart, J. R Gaebler, and D. S. Jin, Nature (London) 454, 
744 (2008). 

[75] J. R Gaebler, J. T. Stewart, T. E. Drake, D. S. Jin, A. Perali, 
R Pieri, and G. C. Strinati, Nat. Phys. 6, 569 (2010). 

[76] L. R Kadanoff and R C. Martin, Phys. Rev. 124, 670 (1961). 

[77] B. R. Patton, Phys. Rev. Lett. 27, 1273 (1971). 

[78] Q. J. Chen, Generalization of BCS theory to short coherence 
length superconductors: A BCS-Bose-Einstein crossover sce- 
nario, Ph.D. thesis, University of Chicago (2000), (freely ac- 
cessible in the ProQuest Dissertations & Theses Database on- 
line). 

[79] Y. He, C.-C. Chien, Q. J. Chen, and K. Levin, Phys. Rev. B 76, 
224516 (2007). 

[80] The contribution of the particle-hole channel T-matrix appears 
(as the last diagram) in Fig. 2.2(a) of Ref. 78, but was later 
dropped in the derivation of the pairing fluctuation theory, 
which focuses primarily on the pairing phenomena in particle- 



particle channel. Nevertheless, one can see from this figure that 
the mixing of Go and G appears in both the particle-particle 
and particle-hole susceptibilities. 
[81] H. Hu, R D. Drummond, and X. J. Liu, Nat. Phys. 3, 469 

(2007) . 

[82] H. Hu, X.-J. Liu, and P. Drummond, Phys. Rev. A 77, 

061605(R) (2008). 
[83] R. Haussmann, W. Rantner, S. Cerrito, and W. Zwerger, Phys. 

Rev. A 75, 023610 (2007). 
[84] S. Floerchinger, M. Scherer, S. Diehl, and C. Wetterich, Phys. 

Rev. B 78, 174528 (2008). 
[85] K. B. Gubbels and H. T. C. Stoof, Phys. Rev. Lett. 100, 140407 

(2008) . 

[86] E. Burovski, N. Prokof'ev, B. Svistunov, and M. Troyer, Phys. 
Rev. Lett. 96, 160402 (2006). 

[87] O. Goulko and M. Wingate, Phys. Rev. A 82, 053621 (2010). 

[88] V. K. Akkineni, D. M. Ceperley, and N. Trivedi, Phys. Rev. B 
76, 165116(2007). 

[89] L. Luo, B. Clancy, J. Joseph, J. Kinast, and J. E. Thomas, Phys. 
Rev. Lett. 98, 080402 (2007). 

[90] L. Luo and J. E. Thomas, J. Low Temp. Phys. 154, 1 (2009). 

[91] M. J. H. Ku, A. T. Sommer, L. W. Cheuk, and M. W. Zwier- 
lein, "Revealing the superfluid lambda transition in the univer- 
sal thermodynamics of a unitary fermi gas," ArXiv: 1 1 10.3309. 

[92] R. Combescot, A. Recati, C. Lobo, and P. Chevy, Phys. Rev. 
Lett. 98, 180402 (2007). 



